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Abstract.  We  present  a large-eddy  simulation  (LES)  of  the  near-field  of 
a low-Mach  number  square  jet.  The  filtered  governing  equations  of  com- 
pressible flows  are  solved  using  the  approximate  Riemann  solver  of  Roe  for 
the  convective  fluxes  and  an  application  of  Gauss’  theorem  for  the  viscous 
fluxes  on  an  unstructured  grid  of  tetrahedral  cells.  The  LES  is  performed 
with  the  Monotone  Integrated  LES  (MILES)  for  the  subgrid  scale  model. 
A square  jet  at  a Mach  number  of  0.3  and  a Reynolds  number  of  3,200  is 
simulated.  Axis-rotation  occurs  and  is  caused  by  Biot-Savart  self-induction 
associated  with  nonuniform  shear-layer  curvature  at  the  exit.  Complex  vor- 
tex topologies  and  interaction  between  vortex  rings  and  rib  vortices  are  ob- 
served downstream.  Mean  properties  are  in  good  agreement  with  the  Direct 
Numerical  Simulation  (DNS)  results  of  Grinstein  et  al.  (1995). 


1.  Introductions 

Turbulent  round  and  plane  jets  are  simple  inhomogeneous  flows  that  can  be 
served  to  verify  models  for  complex  flows  and  have  been  experimentally  and 
numerically  studied  extensively  [1,  2],  Recently,  noncircular  jets  have  been 
gained  much  interest  in  passive  control  due  to  their  enhanced  jet  mixing 
properties  [3,  4,  5,  6,  7]. 

In  the  past,  Reynolds-Averaged  Numerical  Simulations  (RANS)  have 
often  been  employed  for  the  investigation  of  jet  flows.  However,  this  kind  of 
traditional  method  cannot  give  a clear  picture  of  the  unsteady  behaviour 
and  the  results  are  too  sensitive  to  the  turbulence  models  used.  The  advan- 
tage of  Direct  Numerical  Simulation  (DNS)  [8]  and  Large  Eddy  Simulation 
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(LES)  [9]  in  being  able  to  capture  a wide  range  of  spatial  scales  and  tem- 
poral scales  is  of  great  value  for  fundamental  and  applied  studies.  LES  is  a 
good  alternative  to  DNS  which  is  more  expensive  for  high  Reynolds  num- 
ber flows.  One  of  the  practical  uses  of  LES  is  to  predict  the  transition  and 
subsequent  turbulent  mixing  (or  spreading)  process  that  occurs  in  spatially 
developing  flows. 

As  is  well  known,  a free  jet  is  subject  to  Kelvin-Helmholtz  instability 
which,  in  turn,  means  that  the  flow  experiences  an  exponentially  grow- 
ing instability  and  results  in  rolling  up  of  vortices,  pairing  and  merging 
processes.  Linear  instability  analyses  have  confirmed  this  growth.  The  sub- 
sequent breakdown  of  the  large-scale  vortical  structures  is  followed  by  the 
advent  of  spreading  by  the  secondary  three-dimensional  instability  mecha- 
nism. Vortex  dynamics  are  also  expected  to  be  important  in  the  control  of 
noncircular  jets  instability,  transition  from  laminar  to  turbulence  and  the 
jet  development  further  downstream.  Grinstein  and  DeVore  [3]  showed  a 
transition  to  turbulence  in  free  square  jets  characterized  by  the  dynamics 
of  vortex  rings  and  braid  vortices  and  observed  a larger  entrainment  rate 
compared  to  round  jets.  The  underlying  mechanism  for  the  enhanced  en- 
trainment is  the  self-induced  axis-rotation  or  axis-switching  resulting  from 
Biot-Savart  deformation  of  vortex  rings  caused  by  non-uniform  azimuthal 
curvature  at  the  jet  exit.  Foss  and  Zaman  [10]  and  Zaman  [11]  analyzed 
the  influence  of  jet  geometry  on  the  characteristics  of  spreading  in  subsonic 
and  supersonic  flows.  They  observed  a substantial  increase  of  jet  spreading 
when  tabs  were  inserted  at  the  nozzle  which  induced  streamwise  vortex 
pairs. 

The  objective  of  this  research  is  to  validate  the  present  MILES  method- 
ology using  an  unstructured  grid  for  a subsonic  square  jet.  First,  the  im- 
plementation techniques  of  solving  the  governing  equations  are  discussed 
briefly  using  an  approximate  Riemann  solver  of  Roe’s  method  and  the  SGS 
model  for  the  unresolved  small  scales  in  the  LES.  Secondly,  a square  jet  at 
a Reynolds  number  of  3,200  and  a Mach  number  of  0.3  is  simulated.  Tem- 
poral evolutions  are  visualized  to  characterize  the  dynamics  of  deforming 
vortex  rings,  ribs  and  their  interactions.  Statistical  quantities  are  quantified 
and  compared  with  the  DNS  results  of  Grinstein  et  al  [3]. 

2.  Governing  equations 

The  present  LES  method  is  based  on  solving  the  3-D  time-dependent  non- 
dimensional  filtered  transport  equations  for  mass,  momentum  and  energy 
and  the  state  equation  of  ideal  gas  for  fully  compressible  flows  in  an  unstruc- 
tured grid  of  tetrahedral  cells  using  a cell-centered  finite  volume  formation. 
Following  are  the  nondimensionalized  governing  equations  and  the  reference 
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quantities  are:  velocity  U0 0,  static  temperature  T^,  density  Poo,  length  scale 
L and  viscosity  p^,  resulting  in  a Mach  number  of  = Uoo/ and 
a Reynolds  number  of  Re  = Poo^oo-b/Voo- 


d(puj)  djpUjUj)  _ _dp_  + drjj 

dt  dxj  dx,  dxj 

djpe)  d (pe  + p)  Uj  d / 
dt  dxj  dxj  v 


dp  d . 

at  + aTi^)  = 0 

mi  (>  = 1,2,3) 


+ 
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+ qj  + (Tij  + (Tij)Ui) 
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(4) 


where  summation  over  the  three  coordinate  directions  is  implied  in  terms 
with  repeated  indices,  u,  = pu^/p,  and  e = pe/p  are  the  Favre-filtered 
velocity  and  total  energy,  p is  the  pressure,  pe  = p/(7— \)+\u(.Uk  is  the  total 

energy,  qj  = RePr^_X)M a is  the  molecular  heat  flux,  where  Pr  = 0.72 

is  the  molecular  Prandtl  number.  The  dimensionless  molecular  viscosity 
is  dependent  on  temperature,  p = Tn,  where  n = 0.76.  The  subgrid jscale 
(SGS)  stress  is  = pUiUj-pUi  Uj,  the  SGS  heat  transfer  is  Qj  = p(ujT- 
Uj  T),  and  the  molecular  viscous  stress  is  Oij  = f + (| ^ )]• 

The  SGS  terms  Tij  and  Qj  are  modeled  by  the  MILES  method  (t^  = 0 and 
Qj  = 0)  to  account  for  the  unresolved  small  scale  turbulence  [12].  Although 
many  advanced  SGS  models  have  been  developed,  the  simple  models  such 
as  the  MILES  are  widely  used. 

The  inviscid  fluxes  are  computed  by  an  approximate  Riemann  solver 
of  Roe’s  flux  difference  splitting  method  [13],  and  Gauss’  theorem  is  ap- 
plied for  the  viscous  fluxes  and  heat  transfer.  A second-order  Runge-Kutta 
scheme  is  used  for  the  time  marching,  and  the  time  step  is  determined 
such  that  the  total  CFL  (Courant-Friedrichs-Lewy)  number  is  less  than 
1.0.  The  gradients  of  each  variable  (p,  pu,pv,pw,pe)  in  a cell  are  computed 
using  Least  Squares  (LS)  method  of  Singular  Value  Decomposition  (SVD), 
where  eight  neighboring  stencil  cells  are  used.  More  details  of  the  numerical 
procedure  can  be  found  in  Ref.  [14,  15]. 

The  present  MILES  method  has  been  validated  for  several  benchmark 
problems.  In  Okong’o,  Knight  and  Zhou  [14],  good  agreement  was  achieved 
for  the  energy  spectrum  in  an  isotropic  turbulent  flow;  also,  the  comparisons 
with  the  DNS  and  experiment  were  good  for  a channel  flow.  A LES  of  a su- 
personic flat  plate  boundary  layer  was  successfully  validated  by  comparing 
with  experimental  and  theoretical  results  [16]. 
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3.  Details  of  Computation 

A free  square  jet  of  width  D at  Reynolds  number  of  3,200  and  Mach  number 
of  0.3  is  studied.  The  grid  consists  of  65  x 65  x 65  hexahedral  cells  in  an 
unstructured  grid  covering  a computational  domain  of  5 D in  the  streamwise 
direction  and  ±3 D along  the  transverse  directions.  Each  hexahedral  cell  is 
divided  into  five  tetrahedral  cells,  yielding  a total  of  1.3M  tetrahedra.  A 
uniform  grid  is  used  along  the  streamwise  direction  with  the  hexahedral  grid 
spacing  of  Ax/D  — 0.078,  which  is  larger  than  0.04  used  by  Grinstein  et  al. 
in  their  DNS  [3].  The  grids  are  stretched  along  the  other  two  directions  and 
the  minimum  grid  spacings  are  Ay/D  = Az/D  = 0.0375.  The  imposed 
boundary  conditions  include  inflow,  outflow  and  wall  boundaries.  At  the 
inflow,  the  streamwise  velocities  are  prescribed  as 

u = U[  1 4-  A sin(27r/t)]  (5) 


and 


U = 0.5  {Too  [1  - tanh[62(2|j/|/D  - D/(2|y|))]  x 

0.5  UQ o [1  - tanh[M2|z|/D  - D/(2|z|))],  (6) 

where  A — 0.02  is  the  perturbation  amplitude,  / is  the  forcing  frequency 
(/  = 0.5),  62  = O.25iii/0,  where  Ri/0  = 40,  R — D/2  and  0 is  the 
momentum  thickness.  Zero-gradient  condition  is  imposed  at  the  outlet  and 
symmetry  boundary  conditions  are  used  at  the  side  walls. 

4.  LES  Results 

4.1.  VORTEX  DYNAMICS  AND  TOPOLOGY 

The  iso-surfaces  of  the  total  vorticity  u = + u*  + u2z  corresponding 

to  u)  = 0.25wpeafc  are  shown  in  Fig.  1.  Azimuthal  nonuniformities  make 
the  evolution  of  the  jet  shear  layer  more  complicated  relative  to  circular 
jets.  Close  to  the  jet  exit,  a smooth  square  vortex  sheet  can  be  observed, 
and  subsequently  rolled-up  vortex-ring  structures  form  due  to  shear-layer 
Kelvin-Helmholtz  instability.  However,  the  vortex  rings  further  downstream 
deform  to  non-planar  shape  due  to  self-induction  mechanism  caused  by  az- 
imuthal nonuniformities.  The  deformed  vortex  rings  are  connected  with  the 
four  corners  of  the  initial  square  sheet  by  ribs.  The  hairpin  braid  vortices 
aligned  with  the  corners  progress  faster  in  the  diagonal  direction  than  the 
others,  “which  results  in  redistribution  of  energy  between  azimuthal  and 
streamwise  vortices”  [5].  Further  downstream,  the  jet  development  is  char- 
acterized by  the  strong  interaction  between  vortex  rings  and  braid  vortices, 
which  leads  to  a final  breakdown  of  the  large-scales  coherent  structures 
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and  transition  to  the  turbulent  flow.  The  self-induced  deformation  of  the 
rings  and  rib  pair  were  explained  to  be  the  leading  mechanism  for  larger 
entrainment  properties  in  non-circular  jets  relative  to  circular  jets  [7] . “The 
interactions  between  the  streamwise  vortices  and  the  vortex  rings  is  remi- 
niscent of  the  interaction  between  ribs  and  spanwise  rollers  in  the  mixing 
layer”  [7].  Mixing  of  jets  with  surroundings  can  be  enhanced  through  con- 
trolling the  formation,  development  and  interaction  of  large-scale  coherent 
structures  passively. 

Fig.  2 shows  the  instantaneous  crosswise  vorticity  ujz  = dv/dx  — du/dy 
contours  at  a central  x — y plane.  Rolled- up  structures  can  be  observed  near 
the  base,  and  subsequently  symmetrical  counterrotating  toroidal  structures 
form  in  the  shear  layer  which  are  then  followed  by  their  stretching  and 
deformation.  The  organized  structures  can  be  broken  down  into  smaller 
eddies  further  downstream.  Evidently,  the  MILES  model  can  capture  the 
transition  process.  Large-scale  vortex  rings  dominate  in  the  near-field  and 
the  small  vortices  dominate  downstream  after  the  breakdown.  The  spa- 
tial spreading  can  be  clearly  observed  downstream  as  the  jet  spreads  by 
entraining  mass  from  the  surrounding  nonvortical  fluid. 

Fig.  3 (a-d)  shows  the  contours  of  instantaneous  streamwise  vorticity, 
ux  = dv/dz  — dw/dy,  across  the  y — z planes  of  x/D  = 1,  2,  3 and  4. 
Quite  different  behaviour  can  be  observed  at  different  axial  positions.  Vor- 
tex shears  with  some  rounded-corners  are  stretched  and  thickened  but  still 
keep  the  initial  square  shape  at  the  position  of  x/D  = 1.  The  jet  cross  sec- 
tion switches  axis  45°  relative  to  that  of  the  jet  nozzle  at  the  axial  location 
of  x/D  = 2,  3 due  to  self-induced  velocity  around  the  corners  by  the  pres- 
ence of  streamwise  vorticity.  The  flow  structure  develops  into  a irregular 
shape  further  downstream  at  x/D  = 4. 

4.2.  JET  MIXING  AND  ENTRAINMENT 

The  statistical  quantities  are  obtained  by  averaging  over  5 forcing  cycles 
after  a statistically  stationary  state  has  been  reached  after  an  elapsed  di- 
mensionless physical  time  of  10. 

The  centerline  distributions  of  the  mean  axial  velocity  by  LES  compared 
with  the  DNS  results  of  Grinstein  et  al.  [3]  are  plotted  in  Fig.  4.  The  mean 
velocity  initially  decays  within  the  first  1.5  diameters  and  subsequently 
shows  a slight  increase  due  to  the  periodic  roll-up  by  the  sinusoidal  forcing. 
The  decay  after  3.2  diameters  is  the  result  of  turbulent  mixing.  In  general, 
the  agreement  between  present  MILES  and  previous  DNS  results  are  good. 

The  corresponding  centerline  r.m.s.  velocities  u' /Uc  are  shown  in  Fig.  5. 
Note  that  close  to  the  inflow  plane  the  r.m.s.  of  velocity  fluctuations  is 
about  2 percent  and  corresponds  to  the  imposed  disturbance  level.  In  the 
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potential  core  region,  the  r.m.s.  velocity  decreases  slightly  with  increasing 
axial  distance  which  shows  a tendency  to  remain  laminar  with  low  tur- 
bulence intensities.  Beyond  the  end  of  the  potential  core,  the  fluctuations 
of  velocity  increase  very  rapidly,  which  indicates  the  appearance  of  the 
secondary-instability  mechanism  which  leads  to  the  final  breakdown  of  the 
large  vortex  structures.  Due  to  insufficient  length  scale  in  the  streamwise 
direction,  the  self-similar  behaviour  has  not  been  achieved.  However,  it  can 
be  seen  that  the  agreement  between  MILES  and  DNS  is  very  good  for  the 
transient  square  jet. 


5.  Conclusions 

A subsonic  free  square  jet  of  Mach  0.3  and  Reynolds  number  of  3,200  has 
been  investigated  by  the  MILES  using  an  unstructured  grid.  Self-induced 
axis-switching  occurs  due  to  azimuthal  non-uniformity  at  the  exit.  Complex 
vortex  topologies  and  interaction  between  vortex  rings  and  rib  vortices  are 
observed  downstream.  The  mean  properties  are  in  good  agreement  with  the 
DNS  results  of  Grinstein  et  al.  [3]. 
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Figure  1.  Instantaneous  isosurfaces  of  Figure  2.  Instantaneous  streamsise  vor- 
total  vorticity  w = 0.25u>penfc.  ticity  contours  at  the  x — y centre-plane. 


